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Abstract 

The off-diagonal disorder caused by random spin orientations in the para- 
magnetic (PM) state of the double exchange (DE) model is described by 
using the coherent-potential-approximation (CPA), which is combined with 
the variational mean-field approach for the Curie temperature (Tc). Our CPA 
approach is essentially non-local and based on the perturbation theory expan- 
sion for the T-matrix with respect to the fluctuations of hoppings from the 
"mean values" specified by matrix elements of the self energy, so that in the 
first order it becomes identical to the DE theory by de Gennes. The second- 
order effects, considered in the present work, can be viewed as an extension 
of this theory. They are not negligible, and lead to a substantial reduction of 
Tc in the one-orbital case. Even more dramatic changes are expected in the 
case of orbital degeneracy, when each site of the cubic lattice is represented by 
two eg orbitals, which also specify the form of interatomic transfer integrals. 
Particularly, the existence of two Van Hove singularities in the spectrum of 
degenerate DE model (one of which is expected near the Fermi level in the 
30%-doped LaMnOs) may lead to the branching of CPA solutions, when the 
Green function and the self energy become double-valued functions in cer- 
tain region of the complex plane. Such a behavior can be interpreted as an 
intrinsic inhomogeneity of the PM state, which consists of two phases char- 
acterized by different electronic densities. The phase separation occurs below 
certain transition temperature, Tp, and naturally explains the appearance of 
several magnetic transition points, which are frequently seen in manganites. 
We discuss possible implications of our theory to the experimental situation 
in manganites, as well as possible extensions which needs to be done in order 
to clarify its credibility. 

PACS numbers: 75.10.Lp, 75.20.-g, 74.80.-g, 75.30.Vn 
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I. INTRODUCTION 



The nature of paramagnetic (PM) state in perovskite manganese oxides (the manganites) 
is one of the fundamental questions, the answer to which is directly related with understand- 
ing of the phenomenon of colossal magnetoresistance. 

There is no doubts that any theoretical model for manganites should include (at least, as 
one of the main ingredients) the double exchange (DE) physics, which enforces the atomic 
Hund's rule and penalizes the hoppings of polarized Cg electrons to the sites with the op- 
posite direction of the locahzed t2g spins. If the spins are treated classically, the cor- 
responding Hamiltonian is given, in the local coordinate frame specified by the directions 
ei=(cos 0i sin 9i, sin 0i sin 9i, cos ^i) of the spin magnetic moments, by^ 

^ii = (1) 

where are the bare transfer integrals between sites i and j (which can be cither scalars 
or 2x2 matrices depending on how the orbital degeneracy of the Cg states is treated in the 
model), and ^ij=cos ^ cos ^+sin ^ sin ^e~*'^'^'~'^j^ describes their modulations caused by the 
deviation from the ferromagnetic (FM) alignment in the pair i-j. 

The form of Hamiltonian (1) implies that the spin magnetic moments are saturated 
(due to the strong Hund's rule coupling) and the spin disorder corresponding to the PM 
state is in fact an orientational spin disorder. Despite an apparent simplicity of the DE 
Hamiltonian (1), the description of this orientational spin disorder is not an easy task and so 
far there were only few theoretical model, which were based on rather severe (and presumably 
unsatisfactory) approximations. 

The first one was proposed l)y de Gennes more than forty years ago.^ In his theory, all 
^ij are replaced by an averaged value ^, so that the spin disorder enters the model only as 
a renormalization (narrowing) of the e^-bandwidth. The effect is not particularly strong 
and the fully disordered PM state corresponds to ^=|. A generahzation of this theory to 
the case of quantum spins was given by Kubo and Ohata.^ The same idea was exploited 
recently in a number of theories aiming to study the behavior of orbital degrees of freedom 
at elevated temperature,® but based on the same kind of simplifications for the spin disorder 
and its effects on the kinetic energy. 

Another direction, which features more recent activity, is the single-site dynamical mean- 
field theory (DMFT, see Ref. 7 for a review) for the FM Kondo lattice model.^'^ The model 
itself can be viewed as a prototype of the DE Hamiltonian (1) before projecting out the 
minority-spin states in the local coordinate frame.'* If the localized t2g spins are treated 
classically (that is typically the case), this method is similar to the disordered local moment 
approach proposed by Gyorffy et al.,^^ and based on the coherent-potential-approximation 
(CPA) for the electronic structure of the disordered state. 

Now it is almost generally accepted that both approaches are inadequate as they fail 
to explain not even all, but a certain number of observations in manganites of a principal 
character such as the absolute value and the doping dependence of the Curie temperature 
(Tc),^^ the insulating behavior above Tc,^^ and a rich magnetic phase diagram along the 
temperature axis, which typically show a number of magnetic phase transitions^^ and the 
phase coexistence^^"^^ in certain temperature interval. Therefore, it is clear that the theory 
must be revised. 
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Then, there are two possible ways to precede. One is to modernize the model it- 
self by including additional ingredients such as the Jahn- Teller distortion/^ the Coulomb 
correlations,®'^^ and the disorder effects caused by the chemical substitution. Another pos- 
sibility is to stick to the basic concept of the DE physics and try to formulate a more 
advanced theory of the spin disorder described by the Hamiltonian (1), which would go 
beyond the simple scaling theory by de Gennes^ as well as the single-site approximation 
inherent to DMFT.*^'^ 

An attention to the second direction was drawn recently by Varma.^^ The main challenge 
to the theoretical description of the orientational spin disorder in the DE systems comes from 
the fact that it enters the Hamiltonian (1) as an off- diagonal disorder of interatomic transfer 
integrals, which presents a serious and not well investigated problem. In the present work we 
try to investigate some possibilities along this line by employing a non-local CPA approach. 

What do we expect? 

1. It was realized very recently that many aspects of seemingly complicated low- 
temperature behavior of the doped manganites (the rich magnetic phase diagram, 
optical properties, etc.) can be understood from the viewpoint of DE physics, if the 
latter is considered in the combination with the realistic electronic structure for the 
itinerant Cg electrons and takes into account the strong dependence of this electronic 
structure on the magnetic ordering (see, e.g., Ref. 22 and references therein). If this 
scenario is correct and can be extended to the high-temperature regime (that is still a 
big question), there should be something peculiar in the electronic structure of the dis- 
ordered PM state, which can be linked to the unique properties of manganites. In order 
to gain insight into this problem, let us start with the DE picture by de Gennes, and 
consider it in the combination with the correct form of the transfer integrals between 
two Cg orbitals on the cubic lattice. Then, the properties of the PM state should be 
directly related with details of the electronic structure of the FM state, which is shown 
in Fig. 1 and connected with the PM electronic structure by the scaling transformation. 
This electronic structure is indeed very peculiar because of two Van Hove singularities 
at the (tt, tt, 0) and (0, tt, 0) points of the Brillouin zone, which are responsible for two 
kinks of density of states at ±1. It is also interesting to note that the first singularity 
appears near the Fermi surface when the hole concentration is close to 0.3, i.e. in 
the most interesting regime from the viewpoint of colossal magnetoresistance.^^ Such 
a behavior was discussed by Dzero, Gor'kov and Kresin.^^ They also argued that this 
singularity can contribute to the T^/^ dependence of the specific heat in the FM state. 
If so, what is the possible role of these singularities in the case of the spin disorder? 
Note that apart from the single-site approximation, all recent DMFT calculations^'^ 
employed a model semi-circular density of states, and therefore could not address this 
problem. 

2. There are many anticipations in the literature that there is some hidden degree of 
freedom which controls the properties of perovskite manganites. A typical example is 
the picture of orbital disorder proposed in Ref. 19 in connection with the anomalous 
behavior of the optical conductivity in the FM state of perovskite manganites. Ac- 
cording to this picture, the large on-site Coulomb interaction gives rise to the orbital 
polarization at each site of the system. The local orbital polarizations remain even 
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in the cubic FM phase, but without the long-range ordering. In the present work we 
will show that, in principle, by considering non-local effects in the framework of pure 
DE model, one may have an alternative scenario, when there is a certain degree of 
freedom which does control the properties of manganites. However, contrary to the 
orbital polarization, this parameter is essentially non-local and attached to the bond 
of the DE system rather than to the site. 

3. There were many debated about the phase separation in perovskite manganites, ^^'^^ 
and according to some scenarios this effect plays an important role also at ele- 
vated temperatures, being actually the main trigger behind phenomenon of the colos- 
sal magnetoresistance.^° The problem was intensively studied numerically, using the 
Monte Carlo techniques. If this is indeed the case, what does it mean on the lan- 
guage of analytical solutions of the DE model (or its refinements)? Presumably, the 
only possibility to have two (and more) phases at the same time is to admit that the 
self energy (and the Green function) of the model can be a multi-valued function in 
certain region of the complex plane. Such a behavior of non-linear CPA equations 
was considered as one of the main troublemakers in the past,^^'^^ but may have some 
physical explanation in the light of newly proposed ideas of phase separation. 

The rest of the paper is organized as follows. In Sec. II we briefly review the main 
ideas of the variational mean-field approach. In Sec. Ill we describe general ideas of the 
non-local CPA to the problem of orientational spin disorder in the DE systems. In Sec. IV 
we consider the CPA solution for the PM phase of the one-orbital model and evaluate the 
Curie temperature. We will argue that two seemingly different approaches to the problem 
of spin disorder in the DE model, one of which was proposed by de Gcnnes^ and the other 
one is based on the DMFT^'^, have a common basis and can be regarded as CPA- type 
approaches, but supplemented with different types of approximations. In Sec. V we consider 
more realistic two-orbital case for the Cg electrons and argue that it is qualitatively different 
from the behavior of one-orbital model. Particularly, the CPA self-energy becomes the 
double-valued function in certain region of the complex plane, that can be related with an 
intrinsic inhomogcncity of the PM state of the DE model. In Sec. VI we summarize the 
main results of our work, discuss possible connections with the experimental behavior of 
perovskite manganites as well as possible extensions. 



In order to proceed with the finite temperature description of the DE systems we adopt 
the variational mean-field approach. Namely, we assume that the thermal (or orien- 
tational) average of any physical quantity is given in terms of the single spin orientation 
distribution function, which depends only on the angle between the local spin and an effec- 
tive molecular field A, 



while all correlations between different spins are neglected. In the case of PM-FM transi- 
tion, the effective field can be chosen as A=(0, 0, A). 



II. CALCULATION OF THERMAL AVERAGES 




(2) 
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The first complication comes from the fact that the DE Hamihonian (1) is formulated in 
the local coordinate frame, in which the spin quantization axes at different sites are specified 
by the different direction {ei}. Therefore, we should clarify the meaning of orientational 
averaging in the local coordinate frame.^^ Prom our point of view, it is logical that in order 
to calculate the thermal averages associated with an arbitrary chosen site 0, the global 
coordinate frame should be specified by the direction eo=(cos0o sin ^o, sin 0o sin ^^o, cos^o), 
so that at each instant the spin moment at the site is aligned along the ^-direction. The 
averaging over all possible directions eo of the global coordinate frame in the molecular field 
A can be performed as the second step. 

Then, corresponding distribution function at the site is given by Eq. (2). The distri- 
bution functions at other sites can be constructed as follows. The transformation to the 
coordinate frame associated with the site is given by the matrix: 

(cos^ 0o(cos ^0 ~ 1) + 1 sin0ocos0o(cos^o — 1) — cos0osin0o\ 
sin 00 cos 0o(cos 9o — 1) sin^ 0o(cos 6'o — 1) + 1 — sin 0o sin Oq 
cos 00 sin ^0 sin 00 sin ^0 cos^o / 

In the new coordinates, the i-th moment 

has the direction e.=/?ei=(cos 0- sin^^(, sin0. sin cos^Q, and the effective field becomes 
X'=R\={— cos 00 sin ^^o, — sin 0o sin ^^o, cos 9o)X. 

Obviously, while A' and e- depend on ^o and 0o, the new distribution function 
Pi(e;)ocexp(A' • e-) - does not, and up to this stage the transformation to the new coor- 
dinate frame was only the change of the notations. In order to obtain the total distribution 
function -Pi(e.) for i^O, formulated in the local coordinates of the site and taking into 
account the motion ofeo in the molecular field A, Pi{e[) should be averaged over eo with the 
weight po(eo): 

Pi(e;, A) = ^ y dno exp (A' • e'; + A • eo) . (3) 
The normalization constant u is obtained from the condition: 

J dn[P,{e[, X) = 1. 

The form of Eq. (3) implies that the directions of magnetic moments are not correlated (in the 
spirit of the mean-field approach) and the averaging over eo can be performed independently 
for different sites of the system. For the analysis the PM state and the magnetic transition 
temperature, it is sufficient to consider the small-A hmit. Then, Eq. (3) becomes: 

Pi(e;,A)^i^(l + ^cos^;A^). (4) 

The thermal average of the function F{e[), formulated in the coordinate frame of the site 
and taking into account the motion of both eo and ei, is given by 

F(A)= |dO;Pi(e;,A)F(eO. 

The spin entropy can be computed in terms of the molecular field A as:^ 

-TS{X) = ksT J dO'jPi(e;, A)lnPi(e'i, A). 
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In the second order of A this yields (both for i=0 and i^O): 

-TS{X) (5) 

Then, the free energy of the DE model is given by^'^*^ 

T{T,X)^En{T,X)-TS{X), (6) 
where Ed{T, A) is the electron free energy (or the double exchange energy): 

/+00 
dzfT{z-ii)n{z,X), (7) 
-oo 



calculated in terms of the (orientationally averaged) integrated density of states n{z,X). 
/t(2;— //.) = [exp(|^) + l]^^ is the Fermi-Dirac function (/x being the chemical potential). 

The best approximation for the molecular field A is that which minimizes the free energy 
(6). Assuming that the transition to the FM state is continuous (of the second order), Tc 
can be found from equation: 



d'J^(Tc,X) 



dX^ 



= 0. (8) 

A=0 



In practice, the derivative d^EniT, X)/dX^ near A=0 can be calculated using the variational 
properties of n{z, X) in CPA and the Lloyd formula.^^'^^ 



III. NON-LOCAL CPA FOR THE DOUBLE EXCHANGE MODEL 

In this section we discuss some general aspects of the non-local CPA to the problem of 
orientational spin disorder in the DE model specified by the Hamiltonian (1). We attempt to 
describe the disordered system in an average by introducing an effective energy-dependent 
Hamiltonian, 

Hij(^) = E„(^)<5ij-Ey(^)(l-<5ij), (9) 

where Ey is the non-local part of the self energy, which is restricted by the nearest neighbors; 
and Eii is the local (site-diagonal) part. The non-local formulation of CPA is essential 
because in the low-temperature limit TYij should be replaced by the conventional kinetic 
term in which Eij plays a role of the bare transfer integral ty. Therefore, Eij cannot be 
omitted. On the other hand, E^ will be needed in order to formulate a closed system of 
CPA equations. 

Wc require the effective Hamiltonian (9) to preserve the cubic symmetry of the sys- 
tem and be translationally invariant. The first requirement has a different form depending 
on the degeneracy of the problem and the symmetry properties of basis orbitals, and will 
be considered separately for the one-orbital and degenerate DE models. In any case, us- 
ing the symmetry properties, all matrix elements of the self energy {En, Ey} on the cubic 
lattice can be expressed through {Eoo,Eoi} for one of the dimers (for example, 0-1 in 
Fig. 2). Then, the Hamiltonian (9) can be Fourier transformed to the reciprocal space. 
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'Hq{z)—J2j^ the first equation for the orientationally averaged Green 

function can be written as 

^ij(^) = 7^ / ciqe^'^-^^'-^^) [z - nM , (10) 

where the integration goes over the first Brillouin zone with the volume JIbz- 

Unfortunately, the non-local form of the Hamiltonian (9) in the combination with the 
translational invariance do not necessarily guaranty the fulfillment of causality principle, 
that is the single-particle Green function should be analytic in the upper half of complex 
energy plane and satisfy a certain number of physical requirements. This is still 
a largely unresolved problem, despite numerous efforts over decades. We do not have a 
general solution to it either. What we try to do here is simply to investigate the behavior 
of this particular model and try to answer the question whether it is physical or not. Note 
also that neither dynamical cluster approximation'^'' nor cellular DMFT method^^, for which 
the causality can be rigorously proven, can be easily applied to the problem of off-diagonal 
disorder. 

The calculations near Van Hove singularities in the case of degenerate DE model requires 
very accurate integration in the reciprocal space. In the present work we used the mesh 
consisting of 374660 nonequivalent q-points and corresponding to 258x258x258 divisions of 
the reciprocal lattice vectors. 

In order to formulate the CPA equations we consider only site-diagonal and nearest- 
neighbor elements of Gy. Again, using the symmetry properties they can be expressed 
through {Goo, Goi}- In addition, there is a simple relation between Goo and Goi for given 
Soo and Sqi: 

Goo{z) [z - Eoo{z)] + 5^Goi(-2)Sio(z) = 1, (11) 

i 

which follows from the definition of the Green function (10) and the Hamiltonian (9).^^ 
Using this identity, some matrix elements of the Green function can be easily excluded from 
CPA equations. 

In order to obtain the closed system of CPA equations which connects {Sqo, Sqi} with 
{Goo, Goi} we construct the T-matrix:^^ 



n - n{z) 



[i-G(z) ]^n-n(z)J^ \ (12) 



and require the average of scattering due to the fluctuations ATi—Tl—Ti to vanish on every 
site and every bond of the system, i.e.:^^'^"^ 

Too(z)^Toi(z)^0. (13) 

The hat-symbols in Eq. (12) means that all the quantities are infinite matrices in the real 
space and the matrix multiplications imply also summation over the intermediate sites. Note 
that our approach is different from the so-called cluster-CPA^^'^^'^^'^'^ because the matrix 
operations in Eq. (12) are not confined within a finite cluster (the dimer, in our case). 
We believe that our approach is more logical and more consistent with the requirements 
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of the cubic symmetry and the translational invariance of the system, because all dimers 
are equivalent and should equally contribute to the averaged T-matrix. This equivalence is 
artificially broken in the cluster-CPA approach, which takes into account the contributions 
of only those atoms which are confined within the cluster. However, our approach also causes 
some additional difficulties, because the non-local fiuctuations tend to couple an infinite 
number of sites in Eq. (12). Therefore, for the practical purposes we restrict ourselves by 
the perturbation theory expansion up to the second order of A7i: 



f{z) 



n - n{z) 



n - n{z) 



G{z) 



n - n{z) 



(14) 



As we will show, the first term in this expansion corresponds to the approximation considered 
by de Gennes,^ and the next term is the first correction to this approximation. In order to 
evaluate the matrix elements Too{z) and To\{z) in the approximation given by Eq. (14) it 
is necessary to consider the interactions confined within the twelve-atom cluster shown in 
Fig. 2 (obviously, an additional term in the perturbation theory expansion for the T-matrix 
would require a bigger cluster). All such contributions are hsted in Table I. 



IV. ONE-ORBITAL DOUBLE EXCHANGE MODEL 

In the one-orbital case, the effective DE Hamiltonian takes the following form, in the 
reciprocal space: 

7Yq(z) = Soo(^) - 2(c^ + Cy + 02)^01(2;), 

where T^qo{z) and ^01(2;) are C-numbers, c^=cosg^, and all energies throughout in this 
section are in units of the effective transfer integral to, which is related with the Cg-bandwidth 
W in the FM state as to—W/12. Elements of the Green function, Goo{z) and ^01(2;), are 
obtained from Eq. (10). In subsequent derivations we will retain both Goo{z) and G'oi(-z), but 
only for the sake of convenience of the notations, because formally Goi(^) can be expressed 
through Goo{z) using identity (11). 

The self-consistent CPA equations are obtained from the conditions Too{z)=Toi{z)—0. 
In the local coordinate frame associated with the site 0, eo=(0, 0, 1) and all contributions 
to Too{z) and Toi(-2) shown in Table I should be averaged over the directions of magnetic 
moments of remaining sites of the cluster with probability functions given by Eq. (4). This is 
a tedious, but rather straightforward procedure. In order to calculate the thermal averages, 
it is useful to remember identities hsted in Ref. 41. We drop here all details and, just for 
reader's convenience, list in Appendix A the averaged values for all contributions shown in 
Table I without the derivation, so that every step can be easily checked. We also introduce 
short notations for the self energies: Eoo='7o ^-nd Soi=2/3+o"i, and for the Green function: 
Goo=go and Goi=gi- Then, the CPA equations Tqo{z)=0 and Toi{z)=0 can be written in 
the form (for a=0 and 1, respectively): 

(7« = $a((7,^) + *a((7,^)A^ (15) 

where ^ 

^0(0-, fif) = ( ctq + 6ai -) 5fo - 12(To(Ti5fi, 
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/ 13\ 
g) = 2(70(71^0 - (^(^o + ^^'^i + 2(^1 + 



9i, 



and 

T / A 2 4 /22 7 \ 

These equations should be solved self-consistently in combination with the definition (10) 
for the Green function. Different elements of Gi^{z) and ^i}{z) obtained in such a manner 
for the PM state (A=0) are shown in Fig. 3. We note the following: 

1. both Gij{z) and ^ij{z) are analytic in the upper half of the complex plane; 

2. ImGoo{z)<0 in the upper half plane; 

3. the (numerically obtained) integrated density of states lies in the interval 0<n{n)<l 
and takes all values within this interval as the function of chemical potential /i, meaning 
that our system is well defined for all physical values of the electronic density. 

We believe that for our purposes, the fulfillment of these three causality principles is quite 
sufficient and some additional requirements do not necessary apply here.^^ Note also that in 
the one-orbital case there is only one CPA solution in the complex energy plane. 

In first order expansion for the T-matrix with respect to ATi (thereafter all quantities 
corresponding to such approximation will be denoted by tilde), we naturally reproduce 
parameters of the DE model by Gennes:^ $o=$i=^o=0 and \l/i=2/45. The corresponding 
elements of the Green function are also shown in Fig. 3. Although the second-order approach 
gives rise to the large matrix elements of the self energy, the elements of the Green function 
obtained in the first and second order with respect to ATi are surprisingly close (apart from 
a broadening in the second-order approach, caused by the imaginary part of the self energy), 
meaning that there is a good deal of cancellations of different contributions to Gqo and Gqi- 

However, it is not true for the Curie temperatures. While in the first order Tq is solely de- 
termined by Goo, in the second-order it exphcitly depends on both Gqo and matrix elements 
of the self energy, which make significant difference from the canonical behavior.^ 

Indeed, Tq can be obtained from Eq. (8). In order to evaluate the DE energy, we 
start with the PM solution (A=0) and include all contributions of the first order of 
as a perturbation. Employing variational properties of the averaged integrated density of 
states,^^'^^ n{z, A) can be found using the Lloyd formula: 

n{z, A) ~ n{z, 0) + -Im {^o{z)go{z) - 6^i{z)gi{z)} \\ (16) 

TT 

where ^'o, ^l/i, go, and gi correspond to the PM state. '^"^ Moreover, gi{z) can be ex- 
pressed through go{z) using identity (11). Then, the DE energy takes the form: 
Ed{T, X)c^Ed{T, 0)+P(r)A2, where 
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Taking into account the explicit expression for the entropy term, Eq. (5), we arrive at the 
following equation: 

keTc = -QV{Tc). (18) 

In the first order with respect to A7Y (corresponding to the choice o'o=cri=\^o=0 and 
\l/i=2/45) and after replacing /^(z— /x) by 8(2;— /i), ViT) can be expressed through the DE 
energy of the PM state: V^-^Eoi^M). Taking into account that ^i5(PM)=|Ei5(FM), 
we arrive at the well known expression obtained by de Gennes:^ /csTc^— ^£'i)(FM), where 
£'d(FM) is the DE energy of the fully polarized FM state. 

The results of calculations are shown in Fig. 4.'^^ A more accurate treatment of A7Y 
in the expression for the T-matrix significantly reduces Tc (up to 20% in the second-order 
approach). The values Tc~0.20io and 0.17to obtained correspondingly at n=0.5 and n=0.25 
are substantially reduced in comparison with the (local) DMFT approach.^^ We also note a 
good agreement with the results by Alonso et al.,^^ who used similar variational mean-filed 
approach supplemented with the moments-method for computing the averaged density of 
states in the one-orbital DE model. Using the value W^AeV obtained in bands structure 
calculations for the FM state, Tc can be roughly estimated as Tc<800K. The upper bound, 
corresponding to n—0.5, exceeds the experimental data by factor two. Tc can be further 
reduced by taking into consideration the antiferromagnctic (AFM) superexchange (SE) in- 
teractions between the localized spins'^'^^'"^^'^^ and spatial spin correlations. For example, 
according to recent Monte Carlo simulations, the latter can reduce Tc up to 0.12to~460K 
at n=0.5.^^ 



V. DEGENERATE DOUBLE EXCHANGE MODEL FOR THE Cg ELECTRONS 

In this section we consider a more realistic example of the double exchange model in- 
volving two Cg orbitals, which have the following order: |l) = |x^— y^) and |2) = |3z^— r^). 
Then all quantities, such as the transfer integrals tij, the averaged Green function Gij(z), 
and the self energy ^ij{z) become 2x2 matrices in the basis of these orbitals (thereafter, 
the bold symbols will be reserved for such matrix notations). The DE model is formulated 
in the same way as in the one-orbital case. Namely, modulations of the transfer integrals 
are described by Eq. (1) with the complex multipliers .^jj. What important is the peculiar 
form of the t;j matrices on the cubic lattice, given in terms of Slater-Koster integrals of the 
dda type.^^ For example, for the 0-1 and 0-2 bonds parallel to the 2;-axis (see Fig. 2) these 
matrices have the form: 

toi = to2=(o J). (19) 

In the other words, the hoppings along the ^-direction are allowed only between the 
Sz'^—r'^ orbitals. Throughout in this section, the absolute values of the parameter 
(dda) =W/ 6^0. 7 eV will be used as the energy unit. 

The remaining matrix elements in the xy-plane can be obtained by 90° rotations of the 
0-1 bond around the x- and y-axes. Corresponding transformations of the Cg orbitals are 
given by: 
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1 -v^ 
V3 -1 



) 



(20) 



and 




1 ^3 
-1 



(21) 




.23 



(22) 



(23) 



for the hoppings parallel to the y-axis. 

We remind these symmetry properties as an introduction to the analysis of the self energy 
in the case of the spin disorder, which obeys very similar symmetry constraints. Namely, 
since the cubic symmetry is not destroyed by the disorder, the x^—y^ and Sz^— orbitals 
belong to the same representation of the point symmetry group. Therefore, the site-diagonal 
part of the self energy (as well as of the Green function) will be both diagonal and degenerate 
with respect to the orbital indices, i.e.: 



where Sqq=Sqq. On the other hand, the matrix elements associated with the bond 0- 
1 should transform to themselves according to the tetragonal (6*4^,) symmetry. Since the 
x^—y^ and 3z^— orbitals belong to different representations of the 6*4^ group {a\g and h^g^ 
respectively), the corresponding non-local part of the self energy has the form: 



Note that Eqj(2;) is not necessarily zero. The identity toi=0, which holds for the transfer 
integrals, reflects the hidden symmetry of the ordered FM state. However, there is no reason 
to expect that the same identity will be preserved in the case of the spin disorder. Moreover, 
as we will show below, the condition Y}q-^{z)^^ is indispensable in order to formulate the 
closed system of CPA equations. 

Thus, in the case of the orbital degeneracy, there are three independent matrix elements of 
the self energy: Y^q{z)=Yj^q{z) ^ ^oi('^) ^^"^ ^oi('^)- comparison with the one-orbital case 
we gain an additional non-local parameter, which may control the properties of disordered 
DE systems, even in the case of the strictly imposed cubic symmetry. 

The matrix elements of the self energy in the xy-plane can be obtained using the trans- 
formations (20) and (21), which yield: 
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and 



-"03 



■"05 



Then, the effective Hamiltonian takes the following form, in the reciprocal space: 



Tic 



01 



^01 



)(cy-c^) 



_ 1 /(Sil + 3S22)(c, + c,)+4Silc, 



The orientationally averaged Green function can be then obtained from Eq. (10). Similar 

to the self energy, all matrix elements of the Green function can be expressed through Gqq, 
2 1 22 

Gqi, and Gqi, using the symmetry properties (22) and (23). In addition, they satisfy the 
matrix equation (11), from which one can easily exclude one of the matrix elements. 

The self-consistent CPA equations are obtained from the conditions (13). We drop 
here all details and present only the final result (some intermediate expressions for thermal 
averages of different contributions hsted in Table I are given in Appendix B). We also 
introduce short notations for the self-energy: T,QQ—Ij^Q=ao, Eoi=(7i, and Eo|=2/3-|-(T2; and 

for the Green functions: G'oo=5'0) G'oi=5'i! ^^^^ G^oi=fl'2- Then, the CPA equations can be 
presented in the form (15) with a= 0, 1, and 2 corresponding to the conditions Tqq{z)=0, 
=0, and TqI{z)=0. respectively. Using the second-order expression for the T-matrix 



T, 



oil 



Eq. (14), one can obtain the following expressions for the coefficients $q, and 

$0(0-, 9) = (^CTq + 3(7i + 3(72 + ^^9Q- 6(Jo(7i5ri - 6(7o(725'2, 



(24) 



c^o + Y*^? + l^^i + \(^i(^2 + ^(Ji + ^(72 ) gi 



-CTi 



-CTo 



9 

2^i^2 



1 



1 



(25) 



*2(cr,5') = 2(7o(725fo - 



3 3 9 1 1 , 

-di + + -aia2 - -di + -a2 ) gi 

+ + ^^^2 + ^<7l<72 + ^(71 + 2*^2 + ^ ) ^2, 



(26) 



*o((7, 9) = ^0^2 - (y2g^ - -^go, 



(27) 



and 



*i((T, ^) = I ( 7ai + 21(72 + 2) ^1 + (^21ai + 7a2 + ^) 52} , 

2 4 If/ 2\ / 7\ 

^2((^, ^) = ^ ~ 45^°^° ^ 90 I r ''^^ ^ ''^^ ^ 3 j ^ V ^ ^^^^ ^ 3/ 



(28) 



(29) 
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A. CPA solution for the paramagnetic phase 

In this section we consider CPA solutions for the PM phase of the degenerate DE model, 
and argue that the situation is qualitatively different from the one-orbital model, even in the 
case of cubic symmetry. 

Let us first consider the limit |z|^oo. In this case, the matrix elements of the Green 
function have the following asymptotic behavior:'^^ go^z)^^, gi{z)^^, and g2iz)—^-^; and 
the CPA solution in the second order of ^ can be easily obtained analytically from Eqs. (24)- 
(26) as (To(2;)— o'i{z)—>-—^, and o'2{z)—>-—-^2- It is a single- valued solution and in this 
sense the situation is similar to the one-orbital case. The analysis is supported by numerical 
calculations for large but finite \z\, and represents a typical behavior when Im(z)>0.75 
(Fig. 5). The Van Hove singularities are smeared due to the large imaginary part of z and 
the self energy (Fig. 6). 

However, when we start to approach the real axis, the situation changes dramatically. 
In the first calculations of such type we fix Im(2;) and solve CPA equations by moving 
along the real axis and each time starting with the self-consistent self energy obtained for 
the previous value of Ke{z). Then, in certain region of the complex plane we obtain two 
different solutions, depending on whether we move in the positive or negative direction of 
Re{z). A typical hysteresis loop corresponding to Im(2;)=0.7 is shown in Figs. 5 and 6. 

The behavior is related with the existence of two Van Hove singularities at the (0, tt, 0) 
and (vr, tt, 0) points of the Brillouin zone, which are responsible for some sort of instability 
in the system. The singularities become increasingly important when z approaches the real 
axis. The positions of these singularities depend on matrix elements of the self energy, and 
given by Re(Eoo— SSqi-I-SqI) and Re(EQo-|-SQ?|^— SSqi), respectively. Therefore, by choosing 
different starting points for the self energy, the singularities can be shifted either 'to the 
right' or 'to the left' (Fig. 6). Since the CPA equations are non-linear, this may stabihze 
two different solutions. 

A more complete picture can be obtained from Fig. 7, where we plot Re(SQj— Sqi) in 
the complex plane. Depending on the location in the complex plane, the CPA equations 
(15) are again converged to either the same or two different solutions. The double- valued 
behavior of RbI'E'^—'EqY) typically occurs within the shaded area. Note that this area is 
the result of numerical calculations which depend on the choice of the starting point. We do 
not exclude here the possibility that our result may be incomplete and that with a better 
choice of the starting conditions the double-valued area may be enlarged. 

Our analysis is limited by lm{z)^0.5. Further attempts to approach the real axis were 
conjugated with serious difficulties: the topology of solutions becomes increasingly compli- 
cated and may include several additional branches, some of which are presumably unphysical. 
At the present stage we do not have a clear strategy of how to deal with this problem. Of 
course, formally the situation can be regarded as the violation of causality principles. ^^'^^ 
Nevertheless, we would hke to beheve that the double- valued behavior of the Green function 
and the self energy obtained for Im(2;)>0.5 does have a physical meaning and is not an 
artifact of the model analysis. That is because of the following reasons: 

1. for both solutions, the Green function satisfies the inequality ImG'QQ<0; 

2. the system is well defined in the whole interval of densities 0<n(/i)<l (corresponding 
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to the doping range 0<a;<l or the position of the chemical potential /x<l) if and only if 
to take into account the superposition of the two solutions. Conversely, by considering 
only one of the solutions (and disregarding the other one as unphysical), the density 
will exhibit a finite jump (a discontinuity) at certain value of fj,, and the system will 
be undefined within the discontinuity range. Presumably, the discontinuity of the 
electronic density would present even more unphysical behavior than the fact of the 
existence of two CPA solutions. 

3. In some sense the appearance of two CPA solutions fits well into the logic of our work. 
As it was discussed before, the non-local part of the self energy in the degenerate 
case is characterized by two different matrix elements and, in comparison with the 
one-orbital case, acquires an additional degree of freedom (which may be related with 
some non-local order parameter attached to the bond of the system). At the same 
time, the positions of Van Hove singularities, which control the behavior of non-linear 
CPA equations, depend on these matrix elements. Therefore, the appearance of several 
CPA solutions seems to be natural. 

Note that Im(2;)^0.5 corresponds to the position of the first Matsubara pole for 
/csTpaO.lSQ. Taking into account the realistic value of the parameter ((i(i(T)~0.7eV,^^ it 
roughly corresponds to TRil200K, which can be regarded as the lowest estimate for the 
temperature for which our analysis is strictly justified. 

Below we discuss possible physical consequences of the existence of two CPA solutions 
in the PM state. Our scenario is based on the following observations (see Fig. 7): 

1. The existence of the branch-point {B in Fig. 7), which forms two physical branches 
of CPA solutions in certain area of the complex plane. The requirement implies that 
there is a continuous path around the branch-point, which connects the points located 
on two different branches. 

2. Both branches of the (multi- valued) Green function and the self energy are analytic 
(perhaps except the branch-point itself and the branch-edges) . The requirement allows 
us to use the standard theorems of the contour integration: for example, the contour 
integral around the branch-point does not depend on the form of the contour, etc. 

Strictly speaking, the second requirement is a postulate which is solely based on results 
of numerical calculations and at the present stage we do not have a general proof for it. 
Assume, it is correct. Then, the physical interpretation of the multi-valued behavior be- 
comes rather straightforward and two CPA solutions can be linked to two PM phases (with 
different densities) corresponding to the same chemical potential /i. The situation has many 
things in common with the phenomenon of inhomogeneous phase separation, which was 
intensively discussed for manganites.^^'^'' The new aspect in our case is that both phases are 
paramagnetic. The position of the branch-point itself can be related with the temperature, 
below which the PM state becomes intrinsicly inhomogeneous. 

B. energy integration and appearance of two paramagnetic phases 

In this section we discuss some aspects of the energy integration in the complex plane, 
related with the existence of two CPA-branches. Let us consider the integral: 
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/ + 00 
dzfT{z-l2)X{z), 
-oo 



(30) 



where X(/i) is a physical quantity, which can be the density of Cg electrons, the double 
exchange energy or the change of either of them, and X{z) has the same topology in the 
complex plane as the self energy shown in Fig. 7. Then, the behavior of the integral (30) 
will depend on the position of the chemical potential /i with respect to the double-valued 
area. Generally, we should consider three possibilities (see Fig. 8 for notations): 

1. //<Ti. In this case there is only one CPA solution and the integral (30) can be 
evaluated by using the standard techniques (see, e.g., Ref. 51 and references therein). 

2. Ti</x<T2. In this case the integral (30) takes two values for each value of the chem- 
ical potential /i: Xi{iJ,), if the integrand is solely confined within one branch; and 
X2{fJ')=Xi{iJ,)+AX{fi), if it is extended to the second branch. The discontinuity 
AX{fi) is given by the contour integral Ci around the branch-point. It is impor- 
tant here that we do not try to define X{z) as a single- valued function by introducing 
the branch-cuts, which by itself is largely arbitrary procedure.^^ Instead, we treat both 
branches on an equal footing, which inevitably leads to the multi-valued behavior of 
X{iJ,). This integral (30) can be replaced by the contour integral C2 spreading in the 
single- valued area plus residues calculated at Matsubara poles Zn—^x+inkBT{2n-\-l): 

X{pi) = [ dzfriz - n)X{z) - 27rikBTj2X{zn). (31) 

In order to calculate Xi and X2, the Matsubara poles should he on the first and second 
branches of X{z), respectively.^^ 

3. //>T2. In this case the integration along the real axis shall be combined with the 
discontinuity AX given by the contour integral Ci. The integration can be replaced 
by Eq. (31). In this case, all Matsubara poles lie on the single branch. Therefore, 
X{fj,) is a single- valued function. 

As an illustration, we show in Fig. 9 the behavior of the averaged density as the function 
of II. For T=0.24, corresponding to Im(2;o)~0.754, which is slightly above the branch-point, 
the first Matsubara pole falls beyond the double-valued area and n{jj,) shows a normal 
behavior when for each value of /i there is only one value of the electronic density. For 
smaller T, n(/i) may take two different values for the same fi. This fact can be interpreted 
as the coexistence of two different phases. Such a behavior typically occurs in the interval 
— 1.4<yU,<— 0.8, which may also depend on the temperature. As it was mentioned before, any 
single phase fails to define the electronic density in the whole interval 0<n(/x)<l because 
of the discontinuity of n(/x). For example, for T=0.23 one of the phases is not defined in 
the interval 0.34<n<0.38, corresponding to the discontinuity at /i=— 1.25, and the other 
one - for 0.46<rT<0.66 corresponding to ;U=— 0.80. The problem can be resolved only by 
considering the combination of these two phases for which n is defined everywhere in the 
interval 0<n(/i)<l. 

T'pfsO.23 (about 1800K), roughly corresponding to the position of the branch-point, can 
be regarded as the transition temperature to the two-phase state (a point of PM phase 
separation) . 
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C. phase coexistence 



In this section we briefly consider the problem of phase coexistence using a semi- 
quantitative theory of non-interactive pseudoalloy. Namely, we assume that the free energy 
of the mixed paramagnetic state is given by 

^mix(z/) = (1 - + yE^S - TS^M, (32) 

where E^j^ and e]^^ are the DE energies of two PM phases (with lower and higher den- 
sities of the Cg electrons, respectively); y is the "alloy concentration"; and S^i^{y) is the 
configurational mixing entropy: 

-TSr^M = keT [y\ny+{l- y) ln(l - y)] . 

Then, the equilibrium concentration, which minimizes the free energy (32) is given by 
|/=[exp(^^^)+l]~^. The difference of the DE energies AE£,=Ej^^ —Ej^^ can be calculated 
using the definition (7) and the formula (31) for the contour integration. Since the contour 
C2 is confined within the single-valued area, AEd is given by the difference of residues at a 
limited number of Matsubara poles: 

AEn = 2mkBT ^ [n^^) {z^) - n^'^ (zn)] . (33) 

Some details of these calculations can be found in Appendix C. The discontinuity of the 
electronic density An—n^'^^ —n^^^ can be calculated in a similar way. 

Results of these calculations are shown in Fig. 10. Particularly, the equilibrium al- 
loy concentration y is close to 0.5, meaning that the difference AE^ is small and the 
main contribution to the free energy comes from the entropy term. The electronic den- 
sity {n) = {l—y)rif^^^ +yn^'^^ averaged simultaneously over the spin orientations and the alloy 
concentrations shows a discontinuity at the edges of the two-phase region, meaning that 
the system is not defined. This behavior is unphysical and caused by the non-interactive 
approach to the problem of phase coexistence. The averaged densities range, for which two 
PM phases may coexist is typically 0.3-0.7 and depends on the temperature. 



D. Curie temperature 

The Curie temperature can be obtained from Eq. (18). In the case of degenerate DE 
model, the function T>{T) is given by:^^ 

V{T) = — Im / dzfriz - //) {m>o{z)g^{z) - 6 [^i{z)g^{z) + ^!2{z)g2{z)]} . (34) 

TT J -00 

Results of these calculations are shown in Fig. 11. For /x< — 1, the magnetic transition 
temperature appears to be lower than the point of PM phase separation (Tp). Taking 
into account that rp~0.23 and using results of Fig. 10, jiK—l roughly corresponds to the 
densities (n)<0.5. In this region, Tq should be calculated independently for two different 
phases. The calculations can be done using Eq. (31). Not surprisingly that different phases 
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are characterized by different Tc's. Therefore, for the degenerate DE model we expect 
the existence of two magnetic transition points. With the coohng down of the sample, the 
transition to the FM state takes place first in one of the phases, characterized by lower density 
(the hole-rich phase). Then, within the interval T^^<T< T(^\ the FM phase continues to 
coexist with the PM phase, persisting in the hole-deficient part of the sample. Taking into 
account that (dd(7)~0.7eV, the difference of two transition temperatures, which depend on 
can be evaluated as 0<T(1^''— T(1^''<650K. Finally, for T<Tq^ the system exhibits the proper 
FM order established in both phases. 

Formally, the opposite scenario when the phase separation occurs below the magnetic 
transition temperature is also possible, and according to Fig. 11 may take place when //>— 1 
((ri)>0.5). However, the quantitative description of this situation is beyond the small-A 
limit, considered in the present work. 

VI. CONCLUDING REMARKS 

We have applied a non-local CPA approach to the problem of orientational spin disorder 
in the double exchange systems, which was supplemented by a mean-field theory for the 
analysis of magnetic transition temperature. Our CPA approach is based on the perturbation 
theory expansion for the T-matrix with respect to fluctuations of hoppings from the mean 
value of the DE Hamiltonian specified by the matrix elements of the (non-local) self energy, 
so that in the lowest (first) order it automatically recovers the main results of the DE 
theory by de Gennes (the bandwidth narrowing in the PM state, expression for the Curie 
temperature, etc.).^ Our main focus was on the correction of this theory by higher-order 
effects with respect to the fluctuations, which were included up to the second order and 
treated in the real space. 

In the one-orbital case, it led to a substantial reduction of Tc (up to 20%). Nevertheless, 
the obtained value of Tq was largely overestimated in comparison with results of Monte Carlo 
calculations,^^ due to limitations inherent to the mean-fleld approach. Therefore, a sensible 
description of spatial spin correlations, beyond the mean-fleld approximation, presents a 
very important direction for the improvement of our model. 

It appeared, however, that even on the level of mean-fleld theories the situation is far 
from being fully investigated, especially if one takes into account the effects of orbital degen- 
eracy and details of realistic electronic structure for the Sg electrons in the cubic perovskite 
structure. Particularly, two Van Hove singularities at the (tt, tt, 0) and (0, tt, 0) points of the 
Brillouin zone may contribute to the properties of DE systems not only in the ordered FM 
state, but also in the case of the spin disorder, above the magnetic transition temperature. 
The singularities lead to the branching of CPA solutions for the PM state, so that for certain 
values of the chemical potential the system consists of two PM phases with two different 
densities. In such a situation, the PM state becomes intrinsicly inhomogeneous, that also 
determines details of the phase transition to the FM state with the decrease of the temper- 
ature. The magnetic transition is characterized by the existence of two transition points, 
and flrst takes place in one of the phases which constitutes only a fraction of the sample. 
The FM transition in the second phase takes place at yet another temperature, which can 
be signiflcantly lower than the flrst one. Both magnetic transitions are continuous (of the 
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second order). However, they occur separately in two different pliases wfiicli already exist 
in the PM region, i.e. above the first magnetic transition point. 

How is this scenario consistent with the experimental behavior of perovskite manganites? 

It is true that there is no clear experimental evidence supporting the existence of an 
individual temperature of the PM phase separation, Tp, and the intrinsic inhomogeneity of 
the PM phase (unless it is caused by external factors such as the chemical and structural 
inhomogeneities and the grain boundaries). In this sense our result can be regarded as 
the prediction. On the other hand, the phase coexistence below Tq is rather common, and 
was observed in a number of experiments. In addition, our result naturally explains 
the appearance of several magnetic transition points in perovskite manganites. Yet, a clear 
difference is that experimentally one (or sometimes both) magnetic transitions are antifer- 
romagnetic (typically either to the A- or CE-type AFM state), and the situation when two 
consecutive transitions go to the FM state is not realized in practice.^'' Presumably, the 
difference is caused by the hmitation of our analysis by the PM and FM states, while in 
reality for (ri)<0.7 (i.e. when several magnetic transition points are expected) the formation 
of the (A-, CE-, and C-type) AFM structures seems to be more natural. From this point 
of view, results of our work present mainly an academical interest: we pointed out at the 
principal possibility of the existence of several magnetic transition temperatures for the DE 
systems, however the type of the ordered state considered at the low temperature was not 
sufficiently general.^^ 

Finally, we would like to discuss briefly some possible extensions of our model. 

1. As it was already mentioned, any realistic description of the phase diagram of per- 
ovskite manganites would be incomplete without the (A-, CE-, and C-type) AFM 
structures. The minimal model which captures the behavior of doped manganites at 
the low temperature is the double exchange combined with the isotropic AFM SE 
interaction between the localized spins (J^). So, the magnetic phase diagram at T=0 
can be understood in terms of the anisotropy of interatomic DE interactions, caused 
by the anisotropy of magnetic ordering and operating in the background of isotropic 
AFM SE interactions.^^ In this picture, the SE interactions were needed to shift the 
reference point simultaneously in all bonds in the direction of the AFM coupling, while 
the variety of the phase diagram itself is described by the anisotropy of the DE interac- 
tions. However, similar shift is not applicable for the magnetic transition temperature, 
or may not be the major effect of the SE interactions. For example, the shift of the 
Curie temperature by J'^, Tc^Tc—2\J^\,^'^ will also affect the mesh of Matsubara 
poles in Eq. (31). It will require the careful analysis of the topology of CPA solutions 
near the real axis, which will present the main obstacle for such calculations. 

2. Our main results regarding the nontrivial topology of the CPA solutions for the PM 

state have been based on the second order of the perturbation theory expansion for 
the T-matrix. Since the effect was so dramatic, it naturally rises the question about 
importance of the higher-order terms. It seems to be a very important problem for 
the future analysis. 

3. The main effect discussed in this work is based on the existence of Van Hove singular- 
ities in the density of states of degenerate DE model. However, the exact position of 



18 



these singularities, or even the fact of their existence in reahstic compounds depend on 
many other factors, such as the Mn(3rf)-0(2p) hybridization, the cation and structural 
disorder, the purity of sample, etc. All these factors may significantly alter conclusions 
of our work, if one try to apply them to realistic compounds. 
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APPENDIX A: ORIENTATIONAL AVERAGE OF THE T-MATRIX: 

ONE-ORBITAL MODEL 

After averaging over the directions of magnetic moments with the distribution function 
given by Eq. (4), different contributions to the T-matrix (see Table I) become: 

ATiooG'oo^'^oo = cTodo, 

^{AHooGoAHio} = ^{AHoiGioAHoo} = -^'Jo'^iQi + ^<JogiX\ 
. . 15 

ATCoi = (7i - -^X!^, 

AHooGoxAHii = (jIqi, 

= = 2 2 

AT^ooG'oo^'^oi = AHo\GxxAHi\ — —croaigo + —aogoX , 



^{ATioiCioAHoi} = 5alg, - g^7l^lA^ 

i 

Y:{AnolG^Ann} = 5 (af + + g^ - ^ (a, - ^) g^X^ 
AHoiG^oAHoi = (al + ^) 9i - (^^^i + I) 9i>^^ 

and 
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APPENDIX B: ORIENTATIONAL AVERAGE OF THE T-MATRIX: 
THE CASE OF TWO ORBITALS 

Due to degeneracy, it is sufficient to consider the contributions to only one site-diagonal 
element of the averaged T- matrix, say Tqq{z). Results of the orientational averaging for 

22 

Tqq{z) will be identical. Then, for different contributions hsted in Table I we have: 



and 



[^{ATiooGoiATiio}]" = [^{ATioiCioATioo}]'' = -3ao (ai^i + as^s) + -ao92X', 



For the bond 0-1, results of the orientational averaging of the diagonal 11 and 22 matrix 
elements will be different. The contributions to Tq-^{z) are given by: 

[A?^r]^^ = (7i, 

[AHooGoiA7^ii]ii = al9^, 



[AHooGooAHoir = [AHoiGiiAHiiY' = -aoa,9o, 

2 

15' 



E{A7^0iGioA7Yoi}]'' = 2(7?^i + 3(71(72^2 - ^CT,92X\ 



and 



E{AHoiGijA7^ji}]i^ = 

ij 

^ I ((Ji 6(Ti(72 9(72 + ^CTl + 2(72^ S'l + 3 ((7^ - 2(7i(72 - ^(Ti + ^a2^ 92^- 

^ I ( 7(71 + 21(72 + 2) ^1 - (7(71 - 7(72 - ^) 52} A^; 
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while similar contributions to Tq-^{z) are given by: 



[AT^ooGoiATiii]^^ = (7o^^2, 



2 



[ATiooGooATioi]'' = [ATioiGiiATiii]^^ = -aoa29o + t^^^o^oA^ 
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[AT^oiGioATioi]'' = ( + ^) ^2 - ^ (4(72 + ^) <72A^ 



and 



E{A7ioiGijAWji}]^^ = 

ij 

^ ja (cTi - 2(71(72 + (72 - ^(7i + ^(72^ Qi + ^9(7^ + 6(7i(72 + CT^ + ^(7i + ^(72^ (/2| + 

1 { (7^1 - 7^2 - ^) ^1 - (7(71 + 1(72 + I) 92} X\ 



APPENDIX C: CHANGE OF THE INTEGRATED DENSITY OF STATES IN 

THE DOUBLE- VALUED REGION 

In this appendix we discuss some practical aspects of calculations of the difference 
n!^'^\z)—n^^\z) between two CPA solutions in the PM states. According to Ducastelle,^^ 
n{z) is given by the following expression: 



n{z) = ^ImTr i In G{z) - In ( T 



n-n(z) 



Giz))\, (CI) 



where Tr is the trace over site and orbital indices, N is the number of atomic sites, and the 
hat-symbols here stand for the matrices in the orbital and atomic coordinates space. 
Let us begin with the first term. In the second order of G^'^\z)—G^^\z) we have: 

In §(2) _ In §(1) (^) ~ 2 1^(2) _ §(1) I 1^(2) ^ ^(1) I ' 

= 2 [p'\z)r - &^'\z)r] {[^«(.)]-^ + [d(^\z)r 

which can be further transformed using the definition of the Green function (10) as 
ln^(2)(^) _ in^(i)(^) ~ 2 1^(2) (^) - H(^)(^)| |2^ - n^^\z) - W(2)(^)| ' . 

The inverse matrix {...}~^ can be calculated in the same way as the Green function (10). 
Then, if Ro{z) is the site-diagonal element of {...}"\ and Ri{z) and R2{z) are the site- 
off-diagonal ones corresponding to the x^—y'^ and ?>z'^—r^ states for the bond 0-1, we can 
write 

-ixr |ln^(2)(^) _ in§(i)(^)| ~ ^Aao{z)Ro{z) - 12 {Aai{z)Ri{z) + Aa2{z)R2{z)} , 
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where Aaaiz)^a(^^z)-a2\z). 

In the second term of Eq. (CI) (the so-called vortex correction) we expand In up to the 

second order of [H — T-C{z)]. This expansion is necessary to preserve the variational properties 
of our CPA formalism,^"^ which is based on the same approximation (14) for the T-matrix. 
Thus, we have 



In 1 



n - n{z) 



G{z) 



n-n{z) 



+ 



u-n{z) 



G{z) 



n - n{z) 



]g{z). 



In order to calculate the thermal average of this expression, we note that 



H-H{z) 



G{z) 



H - H{z) 



n - H{z) 



which immediately follows from the CPA equations (13) under condition (14). Thus, 



In 1 



n - n{z) 



G{z] 



n - n{z) 



G{z), 



and corresponding contribution to the integrated density of states is given by 
1 



N 



-Tr In 1 



n - n{z) 



G{z)) ~(To(^)^o(^) -3{(Ti(^)c/i(2;) + 0-2(^)^2 (^)}. 
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TABLES 



TABLE I. Different contributions to the site-diagonal (Too) and site-off-diagonal elements (Tbi) 
of the T-matrix in the second order with respect to the fluctuations A7^ij=7^ij— T^ij in the real space. 
The column 'comment' is used to explain the position of intermediate sites used in the summation 
in Fig. 2. 



contribution comment 



element Too: 




AHoo 




A 1-/ /""v A nj 
ZA/toO^J^OO^/lOO 






i = 1-6 


Ei{AWoiGioAWoo} 


i= 1-6 


Ei{A7^oiGiiA7^io} 


i= 1-6 


element Tqi: 








AHooGoiAHii 




AHooGoqAHoi 




AHoiGiiAHii 




I] i { AHo i Cio AHoi } 


i = 1-5 


Ei{AHoiGiiAHii} 


i = 7-11 


AHoiG^oAHoi 




Eij{A7^oiGijAWji} 


(ij) = (3,7), (4,8), 1 
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FIGURES 



X 



0.8 0.4 0.0 -0.4 -0.8(;t,;c,o) (0,0,0) {K,n,0) 




0.5 0.4 0.3 0.2 0.1 0.0 (0,JC,0) (jl,7C,3l) (0,71,0) 

DOS 

FIG. 1. Tight-binding density of states and the dispersion of two eg bands in the ferromagnetic 

state (in units of dda transfer integral). The dotted Hne shows the positions of the Fermi level as a 
function of hole-concentration x, which is related with the integrated density of states as x=l—n. 
Note the existence of two Van Hove singularities at (tt, tt, 0) and (0, 7r,0), responsible for the kinks 
of density of states at ±1. The first singularity is located near the Fermi level when x~0.3 (shown 
by arrow). 
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8^ 



10 







6 



FIG. 2. Atomic sites which will contribute to the orientational average of the T-matrix for the 
dimer 0-1, when T{z) is given by the second-order perturbation theory expression - Eq. (14). 
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-4 -2 2 4 -4 -2 2 4 

Re(z) Re(z) 
FIG. 3. Behavior of matrix elements of the Green function and the self energy along the real 

axis in the one-orbital double exchange model (in units of transfer integral to)- The same results, 

but using the first-order expression for the T-matrix with respect to the fluctuations /SH are shown 

by dotted line. 



0.30 
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n 

FIG. 4. Curie temperature Tq (in units of to) of the one-orbital double exchange model obtained 
using the first- and second-order expression for the T-matrix with respect to the fluctuations A7^ 
as a function of averaged electronic density. 
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-2-1012-2-1012 

Re(z) Re(z) 
FIG. 5. Behavior of matrix elements of the self energy in the paramagnetic state of degenerate 

double exchange. For Im(2;)>0.75 there is only one solution, while for Im(2;)<0.75 one can obtain 

two self-consistent CPA solutions in certain interval of Re(2;) by starting the iterations with the 

self energy obtained for the previous value of Ylc(z) and moving either in the positive or negative 

direction of the real axis (shown as a hysteresis) . 
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FIG. 7. Topology of the non-local part of the self energy in the complex plane. The branch-point 
is denoted by B. The projection shows an approximate position of the double-valued area for the 
CPA solutions. 




Y, H Re(z) 
FIG. 8. Energy integration for the degenerate double exchange model, (a): The integral along 

the real axis plus the discontinuity given by the contour integral Ci around the branch-point B. 

(b): An equivalent expression in terms of the contour integral C2 spreading in the single- valued 

area of the complex plane and residues calculated at Matsubara poles. The latter contributions 

are different for two different branches, that is equivalent to the discontinuity term in the scheme 

(a). 
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(A 

FIG. 9. Averaged electronic density as a function of chemical potential for the paramagnetic 
state of degenerate double exchange model. For T<0.23, n{lJ-) may take two different value for the 
same chemical potential which correspond to two different CPA solutions shown in Fig. 5. 




FIG. 10. Pseudo-alloy picture for the two-phase paramagnetic state of the degenerate double 



exchange model: the change of the double exchange energy AE]:)=E^^ - 



-E^\ the equilibrium alloy 



concentration y, and the density of Cg electrons averaged simultaneously over the spin orientations 
and the alloy concentrations {n)={l—y)n^^^ +yn^'^^ as a function of chemical potential /j, for three 
different temperatures (in units of c/dcr-integral) . The superscripts and stand for the phases 
correspondingly with lower and higher electronic densities. 
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0.35 




FIG. 11. Main transition temperatures for the degenerate double exchange model (in units of 
(i(i(T- integral) . Tp is the transition temperature to the two-phase paramagnetic state. The shaded 
area shows the approximate range of the chemical potentials (//) when the paramagnetic state 
becomes intrinsicly inhomogeneous. T^^ and T^^ are the Curie temperatures for two different 
phases (characterized by lower and higher densities of the Cg electrons, respectively). 
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